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Where is the x c0 (2P)? 
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Although the analysis of the BaBar Collaboration prefers J p = + for the X(3915), it is difficult to 
assign the X(3915) to the Xco(2f )■ We show that there is an indication of the Xco(2P) with a mass 
around 3840 MeV and width of about 200 MeV in the Belle and BaBar data for 77 -4 DD. 

PACS numbers: 14.40.Pq, 13.25.Gv 

With the discovery of many new states in the charmonium mass region, one faces the problem of under- 
standing their nature (for a comprehensive review, see Ref. |[l|]). Among these XYZ states, the Z(3930) 
was discovered by the Belle Collaboration in the process 77 — > DD yfl, and later on confirmed by the 
BaBar Collaboration in the same process [3(]. The experimantal results on mass, angular distributions, and 
r(Z(3930) -> 77)£(Z(3930) ->■ DD) are all consistent with the expectation for the Xc2(2P), which is the 
radially excited P-wave tensor charmonium state. The mass was measured to be 3927.2 ±2.6 MeV [4|]. It is 
below the D*D* threshold so that it decays dominantly into the DD pair in a D-wave, and its decay width 
is 24 ± 6 MeV (jj. So far, this is the only unambiguously identified radially excited P-wave charmonium 
state. 

In the same mass region, another structure was reported first by the Belle Collaboration [5I] in the J/ ijjui 
invariant mass distribution in exclusive B — > KcuJ/ip decays, and confirmed later by the BaBar Collabo- 
ration with more statistics |(], 7]. This structure is identified as the same state observed in the same final 
states with similar mass in the 77 — > J/ipoo process M\. It is called X(3915), and its mass and width 
are 3917.5 ± 2.7 MeV and 27 ± 10 MeV, respectively ||4|]. Very recently, a spin-parity analysis has been 
performed for the process X(3915) — > J/ipco by the BaBar Collaboration [9]], and the results suggest that 
the quantum numbers of this state are J p = + . It was therefore identified as the Xco(2P) state J9J] follow- 
ing the suggestion of Ref. [10]. However, assigning the X(3915) as the Xco(2P) state faces the following 
problems: 

lUj], the partial decay width of the X(3915) -> J/ipu 



(1) As pointed out in, for instance, Refs. |1|, 
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would be > 1 MeV if it is produced similarly in f?-meson decays as the well-understood conventional 
charmonium states. This estimate is consistent with the one using r(X(3915) — > 77)23(X(3915) — > 
J/i/jui), which was reported to be (52±10±3) eV, if its spin is 0, by the BaBar Collaboration [90. The 
value is consistent with an earlier measurement by the Belle Collaboration |8fl. Were the X(3915) 
the x c o(2P) state, its width into two photons would be similar to that of the XcO • Given T(x c o — > 
77) = 2.3±0.2 keV, it is reasonable to assume r(x c o(2P) -> 77) ~ 1 keV. Identifying the X(3915) 
with the Xco(2-P) would give the same estimate as before, r(X(3915) — > J/tpuj) > 1 MeV. Such 
a value is at least one-order-of-magnitude larger than the OZI (Okubo-Zweig-Iizuka) suppressed 
hadronic widths of the ijj(2S) and the ^(3770). 

(2) Normally, the dominant decay channels of a scalar meson should be the open-flavor modes, i.e. 
OZI allowed, if the meson lies above the corresponding thresholds. However, the X(3915) was 
not observed in the DD channel, despite the facts that they can couple in an 5-wave and the DD 
threshold, 3730 MeV, is about 200 MeV lower than the X(3915) mass. 2 Furthermore, when there 
is no suppression due to either isospin or SU(3) breaking, from the OZI rule or because of tiny phase 
space, all known hadrons which decay in an 5-wave have widths of order 100 MeV or even more. 
The width of the X(3915), being comparative to that of the D-wave decaying Xc2(2P), seems too 
small. 

(3) As mentioned in Ref. fl, the mass difference between the Xc2(2P) and the X(3915), 9.7±3.7 MeV, 
is too small for the fine splitting of P-wave charmonia. It is one-order-of-magnitude smaller than the 
fine splitting of IP states M Xc2 - M Xc0 = 141.45 ± 0.32 MeV I4fl. Furthermore, it is even smaller 
than the analogous splitting in bottomonium systems, M Xb2 r 2 p) — ^xboi^P) = ^6-2 ± 0.8 MeV [4j] 
even though the Hamiltonian terms responsible for the fine splitting are proportional to l/rriQ, with 
tuq the heavy quark mass. 

So where should the Xco(2P) be? There have been a few lattice calculations on the mass spectrum 
of excited charmonium states. For the Xcoi^P) state, a quenched calculation gives a mass of 4091 ± 



61 MeV 11411 , while it is predicted to be around 4 GeV in full QCD with a pion as large as 1 GeV IU5I1 



However, effects due to light sea quarks / larg e pion masses could be significant for the excited P-wave 
charmonia nearby open-charm thresholds 11611 . Results calculated with a lower pion mass, though still 
around 400 MeV, are presented in Ref. Ill 711 . The authors calculated the mass differences between the 



1 Assuming that the DD modes dominate the decays of the Xc2{2P), T(xc2(2P) —¥ 77) is about 0.2 keV jj]. It is of the same 
order as T(xc2 -> 77) = 0.51 ± 0.04 keV. 

2 There has been a suggestion that the X(3915) couples to the DD ITHl . However, the result for the width obtained from their fit 
to the 77 — > DD data of both the Belle (2l and Babar |Q| Collaborations, 8.1 ± 9.7 MeV, is consistent with 0. This implies that 
the X(3915) only plays a minor role in their fit. 
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excited charmonia and the i] c in order to reduce the systematic uncertainty due to tuning the charm quark 
mass. The results are 972 ± 9 MeV and 1041 ± 12 MeV for the Xco(2P) and the Xc2^P), respectively. If 
the experimental mass of the rj c is used, the mass of the Xco(2-P) would be 3953 ± 9 MeV, and meanwhile 
the Xc2(2-P) would be about 100 MeV heavier than the measured mass. However, one may extract the 
fine splitting between the Xc2(2-P) and the Xco(2P) from their calculation, which is 69 ± 15 MeV. Using 
the experimental mass of the Xc2(2-P), one gets M Xc0 ( 2 p) = 3858 ±15 MeV. We regard this value as the 
most reliable lattice estimate of the mass of the Xco(2-P) obtained so far. Quark model predictions are in 
the same region. For instance, the masses of the Xco(2P) and the x r o(2P) are predicted to be 3916 MeV 



and 3979 MeV in the Godfrey-Isgur relativized quark model II18L 11911. respectively. Shifting the Xc2(2-P) 
to the observed mass and keeping the value of the fine splitting would give 3856 MeV for the mass of the 
Xco(2-P)- A recent quark model calculation using a screened potential 112011 predicts a mass of 3842 MeV for 
the Xco{2P) \ meanwhile, the prediction for the Xc2(2-P), 3937 MeV, agrees with the experimental value. 

In the data of the process 77 — > DD from both Belle and BaBar, there is a broad bump below the narrow 
peak of the Xc2(2-P). The Belle Collaboration fits the DD invariant mass spectrum in the region 3.80 GeV< 
w <4.2 GeV with a Breit-Wigner function for the Xc2(2-P) plus a background function oc w~ a 12J], with 
w the invariant mass of the DD pair. The BaBar Collaboration fits the spectrum starting from the DD 



threshold with a background function oc \fw 2 — mf (w — m t ) a exp[— f3(w — m t )], where the DD threshold 
is represented by m t Q3[] . However, if the Xco(2-P) has a mass around 3850 MeV as estimated above, it 
would lie in the region of the broad bump, and it could get hidden in such fits. In this paper, we try to fit 
both data sets with two Breit-Wigner functions. Our results indicate that the Xco(2-P) could have a mass 
around 3840 MeV. 

We will assume that all the cross sections of the DD production in photon-photon collisions are due 
to resonant structures. This means that we will neglect backgrounds due to non-resonant contribution and 
semi-inclusive DDX with undetected X (X can be soft pions or photons). We use the same Breit-Wigner 
function for the resonances as the one used by the BaBar Collaboration yfl. Taking into account both the 
phase space and Blatt-Weisskopf centrifugal barrier factor Fl, the function reads as 

e ( W \ = (lX L+1 — F ^ w) m 

W w (w 2 -M 2 ) 2 + T 2 (w)M 2, K > 

where L is the orbital angular momentum of the DD pair, p is the momentum at the center-of-mass (cm) 
frame, po is the cm momentum at w = M, which is the Breit-Wigner mass of the resonance, and the 
energy-dependent width is given by 



\ 2L+1 M 

i» = r(^) —fI( w ), (2) 
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FIG. 1: Fit to the BaBar (left) and Belle (right) data separately. The D-mass sidebands have been subtracted from the 
Belle data. The dashed and dotted lines represent the contributions from the Xco(2P) and the Xc2(2P), respectively. 



22|] is F = 1 for an S 



with T being the width of the resonance at rest. The centrifugal barrier factor (21 
wave, and 

y 1 ^J{R 2 p 2 - 3) 2 + 9#V 
for a D wave. The same value 1.5 GeV as used in Ref. [3] will be taken for the "interaction radius" R. 

We fit to the BaBar and the Belle data separately in the region from the DD threshold to 4.2 GeV with 
four parameters: the mass Mq and width To of a + resonance which couples to the DD in an S-wave, and 
two normalization constants Nq and N2 for the scalar meson and the Xc2(2P), respectively. The mass and 
width of the XciQP) are fixed to 3927 MeV and 24 MeV [4J], respectively. There is no interference between 
these two structures because they are in different partial waves. Contrary to the BaBar data, the Belle data 
are not efficiency corrected. Nevertheless, the Belle efficiency only decreases by 10% for an increase of 
the invariant mass^from 3.8 to 4.2 GeV, and there is no fine structure in the efficiency and background 
distributions J2, 



2311 . Furthermore, the D-mass sidebands have been subtracted from the Belle data used in 
our fit. A comparison of the best fit to the data are shown in Fig. Q] 

The fit results are collected in Table HI where the uncertainties only reflect the statistical errors in the 
fit. One sees that the two resonance assumption gives a reasonable fit to both data sets. The large value 
of x 2 /dof for the fit to the BaBar data comes mainly from a few bins where the event numbers are quite 
separated from their neighbors. Comparing the resulting parameters from the two fits, the difference in the 
values of the mass is 2a, and the values of the width and the ratio of the normalization constants are fully 
consistent with each other. The mass is compatible with the lattice estimate for the mass of the Xco(2-P) 
discussed above, and the width is of the right order for an 5-wave strongly decaying hadron. Furthermore, 
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Data sets 


X 2 /dof 


M (MeV) 


r (MeV) 


N 2 /N 


BaBar [3J 


1.63 


3848.1 ±7.7 


229 ± 26 


0.0121 ±0.0034 


Belle [2J 


0.79 


3825.1 ±8.4 


212 ±29 


0.0132 ±0.0046 



TABLE I: Results of fitting to the BaBar and Belle data, respectively. See text for details. 

its main decay channel should be D D since such a structure has never been reported elsewhere. All these 
properties indicate that the broad bump could be assigned as the x cU (2P). The weighted average I4J] of both 
fits gives 

M Xcoi 2P) = (3837.6 ± 11.5) MeV, r Xco(2P) = (221 ± 19) MeV. (4) 

In summary, it is difficult to accommodate the X(3915) to the still missing Xdo(%P)- Being a J p = + 
particle, the X(3915) is probably of exotic nature. There is an indication that the present data of the 
77 — > DD process already contain signals of the Xeo(2-P) with a mass and width around 3840 MeV and 
200 MeV, respectively. More refined analysis of the data with higher statistics is definitely necessary to 
confirm our assertion. In addition to the photon-photon fusion process, searching for the Xco(2_P) in the 



DD distribution can also be done at BES-III using radiative decays of higher vector charmonia 12411 . and at 
the hadron colliders LHC and FAIR. 
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